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Abstract 

The infrared behaviour of the medium modified gluon propagator in non-equilibrium situations 
is studied in the covariant gauge using the Schwinger-Keldysh closed-time path formalism. It is 
shown that the magnetic screening mass is non-zero at the one loop level whenever the initial gluon 
distribution function is non isotropic with the assumption that the distribution function of the gluon 
is not divergent at zero transverse momentum. For isotropic gluon distribution functions, such as 
those describing local equilibrium, the magnetic mass at one loop level is zero which is consistent 
with finite temperature field theory results. Assuming that a reasonable initial gluon distribution 
function can be obtained from a perturbative QCD calculation of minijets, we determine these 
out of equilibrium values for the initial magnetic and Debye screening masses at energy densities 
appropriate to RHIC and LHC. We also compare the magnetic masses obtained here with those 
obtained using finite temperature lattice QCD methods at similar temperatures at RHIC and LHC. 
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Experiments at RHIC (Au-Au collisions at y/s = 200 GeV) and LHC (Pb-Pb collisions at 
y/s = 5.5 TeV) will provide an excellent opportunity to produce a quark-gluon plasma in the 
laboratory. There is no doubt that an energy density larger than ~ 5 GeV/fm 3 will be 
created during these collisions but it is not at all clear that the partons produced following 
the collision will reach equilibrium. The study of the equilibration of the quark-gluon plasma 
is very crucial because it determines the time evolution of all global quantities such as energy 
density, number density etc.. This study also plays a crucial role in determining many of the 
potential signatures for quark-gluon plasma detection at RHIC. The space-time evolution of 
the parton gas for this non-equilibrium situation can be be modeled by solving semi classical 
relativistic transport equations [§, || [|, |^, |6|, f7|, |8[, |9i [TOR. Central to solving the transport 
equations is what goes into the scattering kernels. Perturbative vacuum expressions for 
gluon scattering suffer from severe infrared problems. One loop medium effects in equilibrium 



provide an electric (Debye) screening mass, but not a magnetic screening mass JTT| . Thus one 
cannot use a one loop resummed finite temperature gluon propagator as an approximation 
to the scattering Kernel because of severe infrared problems in the limit that po = and 
\p\ — > 0. To obtain magnetic screening in equilibrium situations requires a non-perturbative 
lattice QCD calculation. What we would like to point out here is that if we use the CTP 
formalism [TJj with an arbitrary non thermal initial Gaussian density matrix then it is 
possible to obtain at one loop a magnetic screening mass as long as the initial gluon single 
particle distribution function f(k x , k y , k z , to) is not isotropic with the assumption that the 
gluon distribution function is not divergent at zero transverse momentum. That is we assume 
that at t — to one can write a Fourier decomposition of the Gluon field in terms of creation 
and annihilation operators. By a Bogoliubov transformation at t = to one can always set 
the pair distribution functions (a\(k,t = t )a^(q,t = t )\) = {a\(k,t = t )a(q,t = t Q )x) = 0. 
Thus the propagator will have the usual vacuum part and a term which depends on the 
initial expectation value of the number density 

2 

£<4(M = t )a{q,t = t )x) = f(k,t )(2n) 3 S(k-q), (1) 

A=l 

and we have summed over the physical transverse polarizations. For / to correspond to a 
physically realizable quantity the number density as well as energy density has to be finite. 
Thus f{k,t) has to go to zero as k — > oo fast enough so that one obtains finite number 
density and energy density. For Gaussian initial value problems, one only needs to know the 



two-point function at t = to- In our following analysis, we will also need to make a quasi- 
adiabatic approximation so that we will assume, for the purpose of determining the initial 
screening masses, that the system is time-translation invariant. Thus the only difference we 
will assume in our Green's functions from the usual thermal ones will be the choice of an 
anisotropic f(k,to) which will replace the usual Bose Einstein distribution function. This 
approximation has been discussed in detail by Thoma and others |14||. In our calculations, 



we will use a simple ansatz for f(k, to) in which the parameters will be chosen to agree with 
known distributions for minijet production at RHIC and LHC Here we are not suggesting 
that this effect replaces the nonperturbative magnetic screening mass, but that the effect we 
are considering is of the same order of magnitude and already cures the infrared problems 
of the transport theory. 

In what follows we will examine the infrared behaviour of the medium modified gluon 
propagator at one loop using the CTP formalism. The purpose of this paper is to study the 
static limit of the longitudinal and transverse self energy of the gluon (Debye and magnetic 
screening masses) and, in particular, to determine, at the one loop level, how the magnetic 
screening mass depends on the initial f(k,to). Although technically the magnetic screening 
mass is defined as the position of the zero of the inverse propagator (i.e. in the limit 
Po — 0; \p\ ~^ m sa'm 2 sc = n(m| c )) P^H , at one loop the limit of the inverse propagator as 
Po — 0, \p\ — > is gauge invariant (independent of £ for general covariant £ gauges), and 
moreover this limit is the one important for controlling the infrared properties of the collision 
kernel in the transport theory. Thus in this paper we will use the second limiting process 
to define the screening masses. At arbitrary momentum the polarization is not in general 
gauge invariant at one loop. To have a gauge invariant approximation at one loop one can 
make a hard momentum loop approximation as discussed in ]14], 18 1. 



In particular we are interested in nonisotropic nonthermal forms for f(k,t ) consistent 
with known minijet production results. Let us consider an expanding system of partons in 
1+1 dimensions. For this purpose we introduce the flow velocity of the medium 

= (coshr/, 0, 0, sinhr/), (2) 

where r\ = |ln |^ is the space-time rapidity and u^v? = 1. We define the four symmetric 



tensors |L5], [I6L ft? 
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Here T^ v is transverse with respect to the flow-velocity but L^ u and D^ u are mixtures 
of space-like and time-like components. These tensors satisfy the following transversality 
properties with respect to p M : 

p^ip) = p^{p) = 0, P,PuC^(p) = 0. (4) 

In terms of this tensor basis the gluon propagator in the covariant gauge is given by: 



G, v {p) = -iT, p {p)G T {p) - iL, v {p)G L {p) - i^(p)G D {p), (5) 

where G T , G L , G D correspond to T, L and D components respectively of the full gluon 
propagator at the one loop level. The last part G D (p) is identical to the vacuum part |15 



and hence we do not consider it any more. There are separate Dyson-Schwinger equations 
for the different components of the CTP matrix Green's functions that do not couple with 
each other. These equations can be written in the form 



G T >\p)} =[G T >\p)} +Y,[G T '\P) 

J J l,k 



Tl T > L (p) it ■ G^(p) 



'iT,L, 



(6) 



Here i,j,k,l = +, — are the CTP contour labels, and suppression of Lorentz and color 
indices in the above equation is understood. 

In the Keldysh rotated representation of the CTP formalism, in terms of retarded, ad- 
vanced and symmetric Green's functions we have instead 



g r,a(p) = G r\a(p) + G R,a(p) ■ n /?U(p) • g r,a(p)- 
The straightforward solution of the above equation is given by: 
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where the self-energy contains the medium effects. Similar but more complicated equations 
are obtained for the resummed symmetric Green's functions 

G T s > L (p) = [1 + 2f(p)} sgnip^G^ip) - G T A \p)\ 



+ (np L ( P ) - (1 + 2f(p))sgn(p )[U T ^(p) - n^(p)]) x G T R ^ L (p) x G T /{p)- 



(9) 



For the purposes of obtaining the correct kernel for the Boltzmann equation, we need 
the medium improved Feynman propagator for the gluon at one loop level which is just one 
component {++} of the matrix Green's function of the CTP formalism. 



Gp{p) = G(p) , can be written as: 



G(p) 



(10) 



= -[G s (p) + G A (p)+G R (p) , 

where where Ga, G r , Gs stands for advanced, retarded and symmetric green's function re- 
spectively. In the above equation the '+' sign refers to the upper branch in the closed-time 
path. Using the relations of the various self energies one finds [|14], [18| : 



G ++ (p) 



P 



ReU R (p) + ymUsip) 



(11) 



(p 2 - ReU R (p)) 2 + (ImU R (p)) 2 ' 
where ReU R (p) and ImU R (p) are the real and imaginary part of the retarded self energy. 
These self-energies have both longitudinal and transverse parts which, in the static limit 
(po — 0, |p| — >• 0), give Debye and magnetic screening masses of the gluon respectively. In 
the above equation Hs{p) is the the symmetric part of the self energy. 

To obtain the infrared behaviour of this propagator we need to find the static limit of the 
gluon self energy for an anisotropic f(k,t) corresponding to the initial distribution function 



expected from the parton model. In a frozen ghost formalism |T^, [19), the gluon self energy 
is obtained from the gluon loop and tadpole loop as shown in Fig. |l| The ghost does not 
contribute to the medium effect in this formalism because the initial density of states are 
chosen to be that of the physical gluons. All the effects of the ghost are present in the 
vacuum. 

The general expressions for the real and imaginary part of the gluon self energy in non- 
equilibrium in covariant gauge for an expanding gluonic medium have been derived in a 

Here we examine the static limit of these self-energies which play 
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crucial roles to obtain a finite collision integral to study equilibration of the quark-gluon 
plasma at RHIC and LHC The general expressions for the real part of the longitudinal and 
transverse self energy of the gluon loop are given by: 
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FIG. 1: One Loop Graphs for the Gluon Self-Energy 



and 
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G(q,p) 
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The general expression for the tadpole loop contribution is given by: 

d 3 q f(q) 
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To simplify these equations in the infrared limit we expand f{q — p) as: f{q — p) = f(q) — 
P " V q f(q) and neglect the higher order gradients. Similarly we expand \p — q\ as: |p — q\ = 
\q\(l — In the static limit (first taking p = then using \p\ — > 0), and in the rest frame 
(wo = 1, w = 0) we obtain from Eq. fll2"|): 



iten^.o(po = 0, |p| -> 0) = -2<7 2 5 afe iV c 
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and from Eq. (16|): 
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Adding both the above equations we get the expression for the Debye screening mass: 



(18) 



(19) 
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(20) 



which is the real part of the longitudinal self energy. Nc = 3 is used. This equation was 
obtained by various authors ||20|| . Similarly in the static limit and in the rest frame we get 
fromEq. ([13]): 



ReUl, R ( Po = 0, \p\ -+ 0) = A fc iV c 



7 



and from Eq. (|17|): 

n?^(Po = 0) = g 2 5 ab N c J . [| _ 1(0 . pn (22) 

Adding both the above equations we get the expression for the magnetic screening mass: 
m 2 g = [ReU T Gl:R ( Po = 0, \p\ -> 0)] + [ReU^ R (p = 0, \p\ - 0)] 



= 39 7 <# [ 1f ' 11 + te ' #) 1 + ^4 ' 11 " ( ' ' #) 11 (23) 

which is the real part of the transverse self energy. For the imaginary part of the gluon 
self energy at one loop we get: ImTlQ L . R = 0, in the static limit. Note that the above 
formula uses the medium part of the self energy containing a gluon loop and a tadpole 
loop which appears in the resummed gluon propagator (Eq. flllD). The expression for the 
Debye screening mass, we obtained (see Eq. (|20|)), is the same as that obtained by various 



authors |20 ] . The expression we obtain here for the magnetic mass for a non-equilibrium 
gluon distribtution function is new (see Eq. ( p3|) ). There is no approximation present in the 
derivation of Eqs. (p0|) and (|23|). The static limit results are gauge invariant. 



For an isotropic gluon distribution function /(|<f|) we get from the eq. (]20 ) 

Qg 2 



m, 



l = -\ dqqf(q) (24) 



where q = \q\ and from eq. (|23| ) 



m 



2 



0. (25) 

Furthermore for the special isotropic case when the system is described by an equilibrium 
Bose-Einstein distribution function for the gluon Eqs. ( p0|) and ( |23| ) give: 

m 2 D = g 2 T 2 , and m 2 g = (26) 



respectively. These results (Eq. (p6[)) are identical to those obtained by using finite temper- 
ature field theory in QCD assuming that the system is in thermal equilibrium [O, |TB|, |55|]. It 



is interesting to note that the magnetic mass is not only zero at one loop level in equilibrium 
(Eq. (p6|) ) but it is also zero for any isotropic non-equilibrium gluon distribution function 
(see Eq. (|25|) ). Only when the distribution function is non- isotropic one gets a non-zero 
contribution to the magnetic screening mass with the assumption that the distribution func- 
tion of the gluon is not divergent at zero transverse momentum (see below) . This result is 
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not particular to QCD but also will be true for QED when the distribution function is non- 



isotropic. This is explicitly calculated in |22| for the QED case where we have shown that 
we exactly get the same formula for the magnetic screening mass in QED as we obtained in 
this paper for QCD (gluon loop) except that N c g 2 oc e 2 . 

Before considering the situation at RHIC and LHC let us consider an example where 
there is momentum anisotropy in the transverse and longitudinal momentum distribution. 
For this purpose we work in the cylindrical coordinate system: (q t , <j), q z ). From Eq. (|20"D we 
get: 



m i = -(#W^<^f^ (2?) 
From this equation we realize that when f(q) is isotropic, the dependence on p drops out 
and we obtain eq. (p4|) . For anistropic / the mass depends on the direction of p. In 
what follows we will assume p is along the transverse direction and give values only for this 
direction. Similar results can be obtained for the longitudinal choice. Assuming p is along 
the transverse direction we find 



= & J dqt Jd<P /4V^^. (28) 



(2tt) 3 7 """J ^ J ' ^ dq t 

Integrating by parts in q t we get: 



,2 



dq t q t J d<pjy^ f{q t AAz) + J d<f> J dq z [\q z \f(q t , q z , <f>)]qt=o]- (29) 



For an equilibrium distribution function of the form f e „ = — , = = = — , = * 

we get from the above equation: 

2 T 2 2 T 2 

™lt = — + — = 9 2 T 2 (30) 

which is the correct result obtained by using finite temperature QCD. 
Similarly, changing to q t , <fi, q z coordinate system we get from Eq. (|23| ) 

m ° = 8^7 dmt J d v ^ [ ^r ■ [1 + ~w~ ] + ^—f^ • [1 - ~w~ ]] - (31) 

When p points in the transverse direction and we again perform partial integration over q t 
to obtain: 



ni 



2 

at 



)J [2 J dq t q t J d<p cos 2 <p J dq j^ ^^ - J d(j) J dq z [\q z \f{q t ,q z ,<f>))qt=o\- (32) 
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For an equilibrium distribution function of the form: f eo = — 1= — == = — , = 1 = 

the above equation gives: 

3g 2 r/ x vr 2 T 2 , 7r 2 T 2 , s 

m '* = (2fp [(47r) -6- " (47r) -6~ ] = °' (33) 

which is consistent with finite temperature QCD results. 

Before proceeding to compute the initial magnetic screening mass at RHIC and LHC 
situations we will adopt a non-isotropic test distribution function to compute the magnetic 
screening mass from the formula given by eq. (]32|) . We choose a non-isotropic test distribu- 
tion function of the form: 

/ = , 1 (34) 

where h is a parameter for non-isotropy. For ft = 1 we get the usual Bose-Einstein distri- 
bution function. Using the above non-isotropic distribution function we plot the magnetic 
screening mass from the eq. ( |52D in Fig. 2. It can be seen from Fig. 2 that for h = 1 
(corresponding to Bose-Einstein distribution function) we get m gt = and for h ^ 1 
(corresponding to non-isotropic distribution function) we get a non-zero magnetic screening 
mass. 

Now we consider the realistic situation at RHIC and LHC. For the situation at RHIC 
and LHC where the parton distribution function at t — to, f(q,to), describes an out of 
equilibrium situation we can compute the value of these screening masses assuming the 
distribution function can be described by the parton model result for minijets. We note 
that to compute the second term in rhs of eqs. ( p9|) and fl32| ) we need to know the form 
of [f(qt, q z )} qt =o and its behaviour at RHIC and LHC. In this paper we are considering the 
minijet distribution function which are computed by using pQCD applicable above qt = 
1(2) GeV at RHIC (LHC) which are obtained by saturation arguments as studied by several 
authors f2M. We mention here that pQCD is not applicable for small q t , for example below 
1(2) GeV at RHIC (LHC). If one calculates the pQCD minijet production the q t distribution 
behaves as: oc q^ a where a ~ 4 for high q t and ~ 2 for low qt. If one applies pQCD at small 
q t the distribution function f(qt,q z ) is singular at q t = 0. However, for very low q t pQCD 
formulas are not applicable and hence it is not obvious that distribution will be singular at 
qt=0 at RHIC and LHC. For qt — the only computation available at the moment is from 



the McLerran-Venugopalan model [37 where it is shown that at qt=0 the gluon distribution 



behaves as a constant with respect to q t and does not behave as q t a . In this cases we may 
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FIG. 2: The transverse component of the magnetic screening mass as obtained from eq. ( p^ ) by 
using a non-isotropic distribution function of the form: / 



as a function of the 



e V"t+ h "^ T - 1 

non-isotropy parameter h. Note that for h = 1 the distribution function becomes Bose-Einstein 



and hence the magnetic screening mass is found to be zero from eq. (32) 



assume that the gluon distribution does not diverge at qt=0 in the realistic situation at 
RHIC and LHC. In such situations where f(q t = 0,y,(f>) is not singular (q z = q t smhy for 
gluon minijets) one can neglect the boundary term: 



[yjQt +Qzf(Qt,Qz, 
-[q t \ sinh y\f(q t ,y 



qt=oo 
Qt=0 



\qt=0 



■[\Qz\f(quqz,(/>)]qt=o 



(35) 
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since for massless minijets 

q z = (ft sinh y = when (ft = for finite y. (36) 

Here y is the momentum rapidity of the minijet parton. This boundary condition is not true 
in general, and in particular not true for a thermal distribution since a thermal distribution 
function: f eq = eQt ^s^jr-i ls divergent at (ft=0. However if the gluon distribution at q t =0 
behaves as a constant at (ft=0 at RHIC and LHC initial situations |37| then our vanishing 



boundary condition should be valid at RHIC and LHC. In any case a non-perturbative anal- 
ysis of gluon distribution at (ft=0 is beyond the scope of this paper. If the gluon distribution 
behaviour at (ft=0 is found to be divergent in any non-perturbaive calculation unlike the case 
in [|37| then the values reported in this paper might change. We have computed the Debye 



and magnetic screening masses in this paper above q t —l(2) GeV at RHIC(LHC) which is 
similar to the calculations done by several authors for the Debye screening mass |3l| where 
they have adopted similar cut-off values for the minijet momentum in their calculations. 

With the above arguments and with the vanishing boundary conditions (eqs. (^) and 
([36])) we get from eq. (p9|) , after changing to the rapidity variables: ^ = dy: 

i 3a, 



m 



J dq t q t J d<p J dy f(q t ,<j>,y), (37) 



where f{qt,4>,y) is the non-isotropic gluon distribution function. For a cylindrically sym- 
metric system we get: 

6a. 



m 



D * ~ -rr 

7T 



J dq t q t J dy f{q t ,y). (38) 



This is exactly the same equation used by several authors [H]] in the context of minijet 
plasma equilibration in heavy-ion collisions at RHIC and LHC. Similarly using the same 
vanishing boundary condition (eqs. Q35|) and (|36]) ) we get for the magnetic screening mass 
from eq. (p2|) : 

for cylindrically symmetric distribution function f(qt,y)- It can be noted that in eqs. ( |38| ) 
and (|39D one should not use equilibrium distribution function or any other distribution 
function which does not obey the vanishing boundary condition as stated in eqs. (|35D and 
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For conditions pertinent to RHIC and LHC we use the minijet gluon distribution function 
to evaluate the Debye and magnetic screening masses. At high energy the minijet cross 
section can be calculated by using perturbative QCD (pQCD). The leading order minijet 
cross section is given by: 

(J jet = j dp t J dy x J dy 2 ^j^J2 x i U/A^pj.) x 2 f j/A (x 2 ,Pt) a i:j ^ k i(s,i,u). (40) 

ijkl 

Here x\ and x 2 are the light-cone momentum fractions carried by the partons i and j from 
the projectile and the target, respectively, / are the bound-nucleon structure functions and 
Hi and y 2 are the rapidities of the scattered partons. The symbols with carets refer to the 
parton-parton cm. system. The frij^ki are the elementary pQCD parton cross sections. 
As we will be considering a gluon system we include the dominant gluon production cross 
sections at the partonic level which are given by: 

(41) 

and 

9a 2 ut us st 
°9™ = -2fl3-p-p-#]- (42) 

Here a s is the strong coupling constant and 

s = X\X 2 s = Ap 2 t cosh 2 ^ — - ^ 2 ^j . (43) 

The rapidities y±, y 2 and the momentum fractions x±, x 2 are related by, 

xi = Pt {e Vl + e y2 )/V7s, x 2 = p t (e^ + e^)/^7s. (44) 

The limits of integrations are given by: 

Pmin <Pt< , ~ HsTs/pt ~ e- yi ) <y 2 < HVs'/pt - e^), (45) 

2 cosh yi 

with 



\yi\ < ln(yG/2p min + \]s/Ap 2 min - 1). (46) 

In the above equations p min is the minimum transverse momentum above which minijet 
production is computed by using pQCD. We multiply the above minijet cross sections by 
a K factor K — 2 to account for the higher order O(a^) contributions. The minimum 
transverse momentum above which the minijets are computed via pQCD are of the order of 
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Pmin ~ 1 GeV at RHIC and ~ 2 GeV at LHC p4 |. These values are energy dependent and 
are obtained from the saturation arguments. We take p m in = 1 GeV at RHIC and 2 GeV at 
LHC for our computations. The minijet cross section (Eq. (|1])) can be related to the total 
number of partons (N) by 

N jet = T(0) a jet , (47) 

where T(0) = 9 A 2 /8ttR\ is the total number of nucleon-nucleon collisions per unit area for 
central collisions ||25|| . Here R A = 1.1 A 1 ^ 3 is the nuclear radius. A rough estimate of the 
initial volume in which these initial partons are formed at RHIC and LHC can be given by: 
Vq = ttR\tq, where the partons are assumed to be spread by a length r = l/p m in- Assuming 
that the partons are uniformly distributed in the coordinate space (but non-isotropic in 
momentum space) we can easily extract a phase-space gluon distribution function of the 



gluon from the total number of gluon minijets from Eq. (f47f ). The initial distribution 
function of the gluon is then given by: 

1 



f(Pt,yi) = ^7W^ dN3et /d S P (48) 

where 



7tR a t 



d 3 p = d 2 p t dp z = p t d 2 p t cosliyi dy x . (49) 

Using the above minijet initial gluon distribution function in Eq. (|38"D and Eq. (|39| ) we 

get: 



m Dt 



T(0) f r r 1 

ft ^a t o J J J s cosn y 1 i - kl 

(50) 

for the Debye screening mass and: 
m lt = 

T(0) r , r , r , 1 



-3Ka s f dp t f dy x f dy 2 - — ^xif i / A (xi,p 2 t )x 2 f j/A (x 2 ,p 2 t )d- ij ^ k i(s,i, 

i J J J scoshyi ttt, 



7r 2 R 2 A r J J J scoshy lijkl 

(51) 

for the magnetic screening mass of the gluon at the one loop level. Note that in the above 
equation a s occurs outside the p t integration and hence a scale has to be defined, at which 
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this coupling constant has to be determined. For this purpose we take a s as a s (< p\ >) 
where the momentum scale < p\ > is defined by: 

1 f j._ ji f _,„. f j„. 2n Pt 

ijkl 



U), 



if f f ^^Pt \ 

<P t >= J d PtPt J d Vl J d yz—~-12 X l fi/A(xi,P 2 t ) X 2 f j/A {x 2 ,p 2 t ) &ij^ k l{s,t,' 

(52) 



where cr je * is defined by the Eq. ([401). 

In this paper we will be using both GRV98 [2(J and CTEQ6M |27j parametrizations for 



the gluon and quark structure functions inside free proton with EKS98 [28|] parametrizations 
for the nuclear modifications. In Fig. ^ we present the results of the initial gluon distribution 
function (see Eq. fl48|)) at RHIC as a function of the transverse momentum of the gluon for 
different values of the rapidities. The rapidity y is related to the longitudinal momentum p z 
via: p z = p t smhy. 

We present the longitudinal momentum distribution of the initial gluon minijet distribu- 
tion function at RHIC in Fig. ^ for different values of p t . 

Using these gluon minijet distribution functions in eqs. (^) and ( |5"TD we get for RHIC: 

m Dt = 116 MeV and m gt = 82 MeV (53) 

and for LHC: 

m Dt = 150 MeV and m gt = 105 MeV (54) 

by using GRV98 structure functions along with EKS98. The coupling constant a s (< p\ >) 
is found to be 0.287 at RHIC and 0.214 at LHC. If higher order contribution to the minijet 
production would have not been taken into account then our results of screening mass would 



have been J K(= 2) times less than the above values. The above masses may be lower bounds 
to the actual values as we have used a lower transverse momentum cutoff for minijets in 
order that pQCD to be applicable. However, the gluon distribution may be dominant at 



lower p t |37j and hence the magnitude of the screening mass may increase if one can include 
the soft partons into the gluon distribution function. The values we reported in this paper 
are for gluon minijet distribution functions at RHIC and LHC with p min greater than 1 and 
2 GeV respectively. 

Let us now look at the equilibrium situation. Note that at one loop level we get (see Eq. 
(E6|)): m 2 D = g 2 T 2 and m 2 = in equilibrium. Since the one loop magnetic mass is zero in 
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1.5 2.5 3.5 4.5 5.5 6.5 7.5 8.5 9.5 10.5 

Pt (GeV) 

FIG. 3: Initial gluon distribution function at RHIC energies as a function of pt 

equilibrium we will compare our results with the magnetic mass which is obtained by using 
non-perturbative methods. The magnetic mass obtained by using non-perturbative methods 
in equilibrium is given by: m 2 g = |(0.255g 2 T) 2 , see [|[ ^9|. Assuming a temperature of about 
500 MeV at RHIC and by using the coupling constant value a s =0.287 at RHIC (see above) 
we get: m D = gT= 950 MeV and m g = V /|(0.255^ 2 T)= 563 MeV. Assuming T = 1000 MeV 
at LHC and using the LHC coupling constant a s = 0.214 we obtain = 1.639 GeV and 
m g = 840 MeV. In obtaining these masses, one has integrated over all momentum ranges of 
the equilibrium distribution functions. For example, if one uses a Bose-Einstein distribution 
function in Eq. (j2lD and then integrates from (p m in then we obtain m;j = 486MeV 

at RHIC for T=500 MeV and a s =0.287. Similarly for LHC one obtains m D = 840 MeV for 
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FIG. 4: Initial gluon distribution function at RHIC energies as a function of p z 



T=l GeV and a s = 0.214. 

Note that these values are of the same order as that obtained by using the non-equilibrium 
distribution functions at RHIC and LHC. Since the gluon distribution function may be 
dominant at lower p t the magnitude of the screening mass might increase if one can include 
the soft partons into the gluon distribution function f|, [10, |3(| [H], [32, 33|, |34]. The values 
we reported in this paper are for gluon minijet distribution functions at RHIC and LHC 
with p m i n greater than 1 and 2 GeV respectively. Note that due to the asymmetry we 
have computed a specific component (p in transverse direction) of the Debye (mot) and 
magnetic (m g .) screening mass. If one computes the values in all directions their values 
may be even higher. Similar situations hold for magnetic screening masses at RHIC and 
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LHC. As the magnetic mass is a non-perturbative calculation at equilibrium and ours is a 
one loop calculation at non-equilibrium, we expect that a non-perturbative non-equilibrium 
calculation might give a higher magnetic screening mass. The argument is similar to the 
study of non-perturbative calculation for Debye screening mass at finite temperature 

To summarize, we have applied the closed-time path formalism to non-equilibrium sit- 
uations in QCD expected at RHIC and LHC energies to study the infrared behaviour of 
the one loop gluon self energy. We have followed a frozen ghost formalism where the initial 
density of states consists of physical gluons and the ghost is only present in the vacuum level. 
In the infrared limit of the gluon self energy we obtain a non-vanishing magnetic screening 
mass of the gluon at one loop level for non-isotropic gluon distribution functions with the 
assumption that the distribution function of the gluon is not divergent at zero transverse 
momentum. At RHIC and LHC we assumed that the gluon distribution is not divergent 
at qt=0 which is supported by the computation done in |[37|| . With this approximation we 
then applied pQCD above qt =1(2) GeV at RHIC(LHC) and obtain a reasonable initial 
non-equilibrium gluon-minijet distribution function. Using this non-isotropic gluon minijet 
distribution function above q t = 1(2) GeV at RHIC(LHC) we predicted the values of the 
magnetic and Debye screening masses at the initial time. 
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